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Abstract 

We show how one can obtain kink solutions of ordinary differential equations with polynomial 
nonlinearities by an efficient factorization procedure directly related to the factorization of their 
nonlinear polynomial part. We focus on reaction-diffusion equations in the travelling frame and 
damped-anharmonic-oscillator equations. We also report an interesting pairing of the kink solu- 
tions, a result obtained by reversing the factorization brackets in the supersymmetric quantum 
mechanical style. In this way, one gets ordinary differential equations with a different polynomial 
nonlinearity possessing kink solutions of different width but propagating at the same velocity as 
the kinks of the original equation. This pairing of kinks could have many applications. We il- 
lustrate the mathematical procedure with several important cases, among which the generalized 
Fisher equation, the FitzHugh-Nagumo equation, and the polymerization fronts of microtubules. 
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I. INTRODUCTION 



Factorization of second-order linear differential equations, such as the Schrodinger equa- 
tion, is a well established method to get solutions in an algebraic manner [lj]. Here we are 
interested in factorizations of ordinary differential equations (ODE) of the type 

u" + ju' + F(u) = , (1) 

where F(u) is a given polynomial in u. If the independent variable is the time then 7 
is a damping constant and we are in the case of nonlinear damped oscillator equations. 
Many examples of this type are collected in the Appendix of a paper of Tuszyhski et al Q ■ 
However, the coefficient 7 can also play the role of the constant velocity of a traveling front 
if the independent variable is a traveling coordinate used to reduce a reaction-diffusion (RD) 
equation to the ordinary differential form as briefly sketched in the following. These RD 
travelling fronts or kinks are important objects in low dimensional nonlinear phenomenology 
describing topologically-switched configurations in many areas of biology, ecology, chemistry 
and physics. 

Consider a scalar RD equation for u(x,t) 

du ^d 2 u . 

a =^ +•*"<«). P) 

where T> is the diffusion constant and s is the strength of the reaction process. Eq (J2J) can 
be rewritten as 

du d 2 u . 

where the coefficients have been eliminated by the rescalings t = st and x = (s/T>) 1//2 x, and 
dropping the tilde. Usually, the scalar RD equation possesses travelling wave solutions u(£) 
with £ = x — vt, propagating at speed v. For this type of solutions the RD equation turns 
into the ODE 

u" + vu + F(u) = , (4) 

where ' = D — The latter equation has the same form as nonlinear damped oscillator 
equations with the velocity playing the role of the friction constant. 

For applications in physical optics and acoustics it is convenient to write the travelling 
coordinate in the form £ = kx — ut = k(x — vt) with kv = uj. This is a simple scaling by k 



2 



of the previous coordinate turning Eq. (jlj) into the form 



u" + V -u> + ±F{u) = 

that can be changed back to the form of Eq. ([!} by redefining J — t and F(u) 
In general, performing the factorization of Eq. (Q) means the following 



(5) 



D-f 2 {u) D-h{u) 



u = . 



(6) 



This leads to the equation 



u" - ^uu' - hu' - f 2 u' + hhu = 
au 

The following groupings of terms are possible related to different factorizations: 

n 

a) Berkovich grouping: In 1992, Berkovich [3] proposed to group the terms as follows 

dfi , 



(7) 



u"-(f 1 + f2)u'+[f 1 f 2 -j±u')u = 0, 



(8) 



and furthermore discussed a theorem according to which any factorization of an ODE of the 
form given in Eq. (jUJ) allows to find a class of solutions that can be obtained from solving the 
first-order differential equation u' = j\u. Substituting the latter expression in the Berkovich 
grouping one gets 

dfu 



U" ~ (fib + f2b) U + (^fuf2b - ~j^f lbU ^j U = 



(9) 



where we redefined fi = fu and f 2 = f2b to distinguish this case from our proposal following 
next. For the specific form of the ODEs we consider here, Berkovich's conditions read 

F(u) 



r I t dfib 

fib [ -7 - fib ~ 



u 



(10) 



fib + hb = -7 • 

b) Grouping of this work. We propose here the different grouping of terms 



u" - 



du 



-U + 01 + 02 ) U + 0102« = 



(11) 



(12) 



that can be considered the result of changing the Berkovich factorization by setting fib = 0i 
and f 2 b — > 2 under the conditions 

. . F(u) 



>1V2 



U 



(13) 



01 + 02 + -3— U = -7 . 

du 

The following simple relationship exists between the factoring functions: 

/ , df lb 

02 = J2b ~~j U 



(14) 



and further (third, and so forth) factorizations can be obtained through linear combinations 
of the functions fib, fib and 2 . 

Based on our experience, we think that the grouping we propose is more convenient than 
that of Berkovich and also of other people employing more difficult procedures. The main 
advantage resides in the fact that whereas in Berkovich's scheme Eq. fllUj) is still a differential 
equation to be solved, in our scheme we make a choice of the factorization functions by 
merely factoring polynomial expressions according to Eq. (jlUjl and then imposing Eq. ()14)) 
leads easily to an n-depending 7 coefficient for which the factorization works. This fact 
makes our approach extremely efficient in finding particular solutions of the kink type as 
one can see in the following. 

We will show next on the explicit case of the generalized Fisher equation all the mathe- 
matical constructions related to the factorization brackets and their supersymmetric quan- 
tum mechanical like reverse factorization. In addition, in less detail, we treat within the 
same approach, damped nonlinear oscillators of Dixon- Tuszyhski-Otwinowski type and the 
FitzHugh-Nagumo equation. 

II. GENERALIZED FISHER EQUATION 

Let us consider the generalized Fisher equation given by 



The case n = 1 refers to the common Fisher equation and it will be shortly discussed as a 
subcase. Eq. (J13J) allows to factorize the polynomial function 




(15) 




(16) 



u 



Now, by choosing 




(17) 
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the explicit forms of a% and 7 can be obtained from Eq. (|14j) 



-p-u + 0i + 2 = --a lM "/ 2 + o x (l - « n/2 ) + (l/ai)(l + « n/2 ) 
Introducing the notation h n = (| + l) 1 / 2 one gets 

ai = ±h~ 1 , 7 = =f (hn + h~ l ) . 

Then Eq. ([To]) becomes 



"7 



" ± (h n + h' 1 ) u' + u(l - u n ) = 



and the corresponding factorization is 

D ± /i„(w n/2 + 1)1 \D T K\u n/2 - 1) 



u = . 



It follows that Eq. (J2"U|) is compatible with the first-order differential equation 

u t K 1 (u n/2 - l) u = . 
Integration of Eq. (j2Z)) gives for 7 > 



uz — 1 1 ± exp 



(^-^(e-eo)]) 



-2/n 



Rewritten in the hyperbolic form, we get 



,11 , 

= tanh 

2 2 



1 1 , 

coth 

2 2 



2/n 



2/n 



(18) 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



The tanh form is precisely the solution obtained long ago by Wang 4J and Hereman and 
Takaoka [5J] by more complicated means. 

Moreover, a different solution is possible for 7 < 

-2/n 



1 ± exp 



or 



1 1 

" ' = 1 2 + 2 tan 



"■ = 1 2 + ^ coth 



- 5 P»» - - Co) 

- ~ P*n - K 1 ) (i ~ 6) 



(25) 



2/n 



2/n 



(26) 



respectively. 

2.1 Reversion of factorization brackets without the change of the scaling factors 
Choosing now <f) X = ai(l + u n ^ 2 ) and 2 = (w n ^ 2 — l) leads to the same equation (|20)l 
but now with the factorization 



D =F h n {u n l 2 — 1) Di/i-^u^ + l) 



it = 



and therefore the compatibility is with the different first-order equation 

v! ± h- 1 (u n/2 + l)u = 0. 



(27) 



(28) 



However, the direct integration gives the solution (for 7 > 0) 



u 



1 



liexpK^-vue-eo)] 



2/n 



'-l) 2/n (1± 



exp 



(^-0(^-^0 



-2/n 



(29) 

which are similar to the known solution Eq. (f2*3*|) . For 7 < 0, solutions of the type given by 
Eq. (|25jl are obtained. 

Direct reversion of factorization brackets 
Let us perform now a direct inversion of the factorization brackets in (}2"T]) similar to 
what is done in supersymmetric quantum mechanics in order to enlarge the class of exactly 
solvable quantum hamiltonians 



D =p h-\u n/2 - 1)1 \D ± /i„(M n/2 + 1) 



u = . 



Doing the product of differential operators the following RD equation is obtained 

u" ± (/i n + h" 1 ) u' + u[l+ u n/2 ] [1 - h A n u n/2 ] = . 

Eq. (|5Tjl is compatible with the equation 

u ± /i n (w n/2 + l) u = 0, 

and integration of the latter gives the kink solution of Eq. (jHlj) 

1 



(30) 



(31) 



(32) 



1 ± exp 



(33) 



liexpf^-^)^-^)], 
for 7 > 0. On the other hand, for 7 < the exponent is the same but of opposite sign. 
Hyperbolic forms of the latter solutions are easy to write down and are similar up to widths 
to Eqs. (|2*lj) and (|2l)|). respectively. 



Thus, a different RD equation given by (|3ip with modified polynomial terms and its 
solution have been found by reverting the factorization terms of Eq. (|20|l. Although the 
reaction polynomial is different the velocity parameter remains the same. This is the main 
result of this work: At the velocity corresponding to the travelling kink of a given RD equation 
there is another propagating kink corresponding to a different RD equation that is related to 
the original one by reverse factorization. We can call this kink as the supersymmetric (susy) 
kink because of the mathematical construction. 

Finally, one can ask if the process of reverse factorization can be continued with Eq. 1)31)). 
It can be shown that this is not the case because Eq. (|31jl has already a discretized 
(polynomial-order-dependent) 7 and this fact prevents further solutions of this type. Sup- 
pose we consider the following factorization functions 



a, 1 



hiu n ' 2 



(34) 



Then, one gets a,\ = ±/i^ and solve + a\ — h^ 1 + h n . The solutions are: n — 0, which 
implies linearity, and n = —4, which leads to a Milne-Pinney equation. On the other hand, 
Eq. (|31|) with an arbitrary 7 can be treated by the inverse factorization procedure to get the 
susy partner RD equation and its susy kink. 
2.3 Subcase n = 1 

This subcase is the original Fisher equation describing the propagation of mutant genes 

du d 2 u 
dt dx 2 



+ u(l - u) 



(35) 



In the travelling frame, the Fisher equation has the form 



u" + av! + u(l - u) = . (36) 

When the 7 parameter takes the value 71 = |-\/6 (i.e., h\ = ^) one can factor Fisher's 
equation and employing our method leads easily to the known kink solution 



Up = - ( 1 — tanh 



L #?-&) 



(37) 



that was first obtained by Ablowitz and Zeppetella |7| with a series solution method. On 
the other hand, the susy kink for this case reads 



susy 



1 — tanh 



V6 



(£-6) 



(38) 



i.e., it has a width one and a half times greater than the common Fisher kink and is a 
solution of the partner equation 

u" + ———u' +u(l- -u 1 ' 2 - -u) = . (39) 



6 V 4 4 

A plot of the kinks up and «F,sus y is displayed in Fig. 1. 
2.4 Subcase n = 6 

This subcase is of interest in the light of experiments on polymerization patterns of mi- 
crotubules in centrifuges. It has been discovered that the polymerization of the tubulin 
dimers proceeds in a kink-switching fashion propagating with a constant velocity within the 
sample. Portet, Tuszynski and Dixon p used RD equations to discuss the modification of 
self-organization patterns of MTs as well as the tubulin polymerization under the influence 
of reduced gravitational fields. They used the value n = 6 for the mean critical number of 
tubulin dimers at which the polymerization process starts and showed that the same nucle- 
ation number enters the polynomial term of the RD process for the number concentration c 
of tubulin dimers 

c" + y + c (1 - c 6 ) = . (40) 



1/3 

(£-£0) ) • ( 4 1) 



The polymerization kink in their work reads 

_i ( r3 
cptd = 2 3 I 1 — tanh - 

On the other hand, the susy polymerization kink (see Fig. (2)) of the form 

1/3 



Csusy 



2^(l-tanh[3(£-£o)]) (42) 



can be taken into account according to the hyperbolic form of Eq. ()33)1 . It propagates with 
the same speed and corresponds to the equation 

c" ± -c' + c (1 - 15c 3 - 16c 6 ) = . (43) 

In principle, this equation could be obtained as a consequence of modifying the kinetics 
steps in the microtubule polymerization process. 



S 



III. EQUATIONS OF THE DIXON-TUSZYNSKI-OTWINOWSKI TYPE 



In the context of damped anharmonic oscillators, Dixon et al |8j studied equations of the 



type (in this section, we use ' = D T = 4-) 

u" + u' + Au-u 11 ' 1 =u" + u +u(VA-u%- 1 )(VA + u%- 1 ) = 



(44) 



and gave solutions for the cases A — | and A = ^, with n = 4 and n = 6, respectively. For 
this case, time is the independent variable. The factorization method works nicely if one 
uses g n = \Jn/2 and dealing with the more general equation 



u" ± VA(g n + g- l )u' + u(A - u n ~ 2 ) = 0, 



for which we can employ either the factorization functions 



(45) 



or 



^g^fVI-u^ x ), 2 = T9n (v 7 ! + u^- 1 



(46) 



(47) 



Then, Eq. (I45j) can be factored in the forms 



d t ± g^- 1 + Va)} \d t t <7" V" 1 - VI) 



and 



D T =F g^- 1 - VA)j [D T ± ff" 1 ^" 1 + v 7 !) 
Thus, Eq. (|43j) is compatible with the equations 



u = 



it = 



/ -i / --i 



A) u = 0, 
u' ± (ut- 1 + v 7 !) u = 



(48) 
(49) 

(50) 
(51) 



that follows from Eq. (}4*%j) and Eq. (fUJj). Integration of Eqs. (J5UJ) . (JBTj) gives the solution of 
Eq. (gSD 

v 7 ^ 



2 

n-2 



1 ± exp 



^(9n -gn X ){ 



7 > 



(52) 



9 



and 



2 

n— 2 



.4 



7 < . 



1 ± exp - - 9 n X ) (T - T ) 

The solutions obtained by Dixon et al are particular cases of the latter formulas. 
Reversing now the factorization brackets in (|48|) 

Dr =F 0n 1 _1 - V 7 !) ± 2n(« f ~* + V 7 !)! M = 



leads to the following equation 
which is compatible with the equation 



2.-1 

«2 



A 1 

4 



whose integration gives the solution of Eq. 



1 ± exp[y/Ag n (r - r )] 



and 



2 

n-2 



2 

n-2 



1 ± exp[-\/Z# n (r - r ) 



7 > 



7 < . 



(53) 



(54) 



(55) 



(56) 



(57) 



(58) 



IV. FITZHUGH-NAGUMO EQUATION 

Let us consider the FitzHugh-Nagumo equation 

du d 2 u 



+ u(l — it) (a — u) = 



(59) 



dt dx 2 

where a is a real constant. If a = —1, one gets the real Newell- Whitehead equation describing 
the dynamical behaviour near the bifurcation point for the Rayleigh-Benard convection of 
binary fluid mixtures. The travelling frame form of (|59jl has been discussed in detail by 
Hereman and Takaoka 



1 1 



u" + 71/ + u(u -l)(a-u) = 0. 



The FitzHugh-Nagumo polynomial function allows the following factorizations: 

0x = ±(V2)- 1 (u-l), 2 = ±V2(a-u) 



(60) 



(61) 
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when the 7 parameter is equal to j a i = ±^±I that we also write as j a> i = ±y / a((7 a i — g~ x ) 
where g a \ = — v2a. 

In addition, we can employ the factorization functions 



>x = ±(v / 2)- 1 ( 



a — u, 



±V2(u - 1) 



(62) 



when 7 aj2 
where g a2 



± ^t 2 , or written again in the more symmetric form 7 a 2 = ±yfa(g a 2 — g^) 
-■s/a/2. Thus, Eq. f)0J) can be factored in the two cases 



and 



u" ± j a iu' + u{u — l)(a — u) = 0, 



u" ± 7 a> 2M' + u{u — l)(a — u) = . 



(63) 



(64) 



In passing, we notice that for the Newell- Whitehead case a = — 1 the two equations coincide 
and are the same as the generalized Fisher equation for n = 2. 

In factorization bracket forms, Eqs. (J63|) and (J64|) are written as follows 



D =f V2(a - u) 



D±(v / 2)- 1 (l-n 



and 



a — u) 

and are compatible with the first order differential equations 



u = 



u = 0, 



u' ± (V2) -1 (l - u)u = 0, for 7a(1 , 
w' =F (v / 2) _1 (a — «)« = 0, for7 a 



7a,2 



(65) 
(66) 

(67) 
(68) 



Integration of the latter equations gives the solution of Eq. (160 jl for the two different values 
of the wave front velocity 7 0i i and 7^2. 
For Eq. (jo^j) we get 

1 1 



liexpKv^)- 1 ^-^)] 
for 7 a l positive and negative, respectively. 
As for Eq. (JMj) . the solutions are 



l±exp[-( v / 2)- 1 (e-6 



(69) 



it- 



l±exp[-( v / 2)- 1 a(e-^o) 



liexp^-Me-Co)] 1 



(70) 
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for 7 a> 2 positive and negative, respectively. 

Considering now the factorizations (J55)l and ([66)1 . the change of order of the factorization 
brackets gives 



D ± (V2)-\l - u) \D =f \/2(a - 



u 



and 



L> =p (V2) _1 (a - w) 



u = 



u = 



(71) 



(72) 



Doing the product of differential operators (and considering the factorization term v! — <p2U 
0) gives the following RD equations 



u" ± 7aiw' + u(Au - I) (a - u) = 0, 



and 



u" ± 7a2w' + u(u — l)(a — u — 3m 2 ) = 0, 
Eqs. (J75J) and (J73J) are compatible with the equations 



u =f V2(a - u)« = 



(73) 



(74) 



(75) 



and 



u =f V2(u - l)u = , (76) 

respectively. Integrations of Eqs. (J75j) and (J75J1 give the solutions of Eqs. (|75j) and (JUj), 
respectively. The explicit forms are the following: 
(i) for (EH)) 



1 ± exp[- v / 2a(e - £„)] ' M< 1 ± exp[ v / 2a(e - Co)] ' 



(77) 



(ii) for (J71D 



1 ± exp[V2(e - Co 



l±exp[-V2(£-£ )] 



(78) 



V. CONCLUSION 



This paper has been concerned with stating an efficient factorization scheme of ordinary 
differential equations with polynomial nonlinearities that leads to an easy finding of analyt- 
ical solutions of the kink type that previously have been obtained by far more cumbersome 
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procedures. The main result is an interesting pairing between equations with different poly- 
nomial nonlinearities, which is obtained by applying the susy quantum mechanical reverse 
factorization. The kinks of the two nonlinear equations are of different widths but they 
propagate at the same velocity, or if we deal with damped polynomial nonlinear oscillators 
the two kink solutions correspond to the same friction coefficient. Several important cases, 
such as the generalized Fisher and the FitzHugh-Nagumo equations, have been shown to be 
simple mathematical exercises for this factorization technique. The physical prediction is 
that for commonly occurring propagating fronts, there are two kink fronts of different widths 
at a given propagating velocity. Moreover, the reverse factorization procedure can be also 
applied to the Berkovich scheme with similar results. It will be interesting to apply the 
approach of this work to the discrete case in which various exact results have been obtained 
in recent years More general cases in which the coefficient 7 is an arbitrary function 
could also be of much interest because of possible applications. 
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